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Abstract. 

Exact and approximate formulas for the upper bound of the relative energy 
difference of two Gaussian states with the fixed fidelity between them are derived. 
The reciprocal formulas for the upper bound of the fidelity for the fixed value of 
the relative energy difference are obtained as well. The bounds appear higher for 
pure states than for mixed ones, and their maximal values correspond to squeezed 
vacuum states. In particular, to guarantee the relative energy difference less than 
10%, for quite arbitrary Gaussian states, the fidelity between them must exceed 
the level 0.998866. 
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1. Introduction 

There are many areas of quantum physics (e.g., quantum teleportation or quantum 
state engineering) where one has to compare somehow different quantum states. In 
particular, in many cases it is important to know, how 'close' are two quantum states 
(e.g., the initial one and its teleported or reconstructed partner), described by the 
statistical operators pi and p2- An accepted measure of closeness is the Bures- 
Uhlmann fidelity [1] 

n 2 



Tr \fhhVh^ 



(1) 



Nowadays one can meet this quantity in almost every paper on quantum information, 
and many interesting results have been obtained for the past decade. For 
example, various boundaries or critical values for fidelities, corresponding to different 
teleportation protocols of some specific classes of quantum states, were derived in 
[5H1]. In particular, the value T — Xjl has been established as the boundary 
between classical and quantum domains in the teleportation of coherent states of 
the electromagnetic field [2j, while another critical value = 2/3 was found in [HH]. 
The meaning of these numbers was further elucidated in [5l|6] . Experimentalists also 
frequently use the fidelity as measure of quality of their achievements. A remarkable 
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progress is observed for the past decade: the reported experimental values increased 
from 0.58 0, 0.62 [S] and 0.64 to 0.87 [TO], 0.93 [H] and 0.98 [H]. 

However, looking at these numbers, the following inevitable question arises: is 
it sufhcient to have, say, 90% fidelity (or even 98%), to be sure that the two states 
are indeed 'similar' or 'close' to each other? It seems that the answer depends on 
the concrete situation and additional information or assumptions about the states, in 
particular, on the exact meaning of the word 'close'. 

The essence of the problem and motivation to study it can be elucidated in the 
following simple example. Consider two coherent states la) and |/3). Their fidelity 
equals — |(a|/3)p = exp (— |a — /3p), so it depends on the difference \a — I3\ only. 
The question is: are the two pairs of states, {|a), |/3)} and {\a + A),\I3 + A)}, 'equally 
close' (or 'equally distant')? From the pure geometrical point of view, the answer is 
definitely positive, since the second pair can be obtained from the first one by means 
of a unitary transformation or by a simple shift in the complex plane of parameters 
a. From this point of view, the space of parameters is homogeneous. But it is 
certainly non-homogeneous from the physical point of view, since the vacuum state |0) 
is obviously distinguished (in this connection, the 'energy-sensitive' distances between 
quantum states were proposed in [13]). Let us take |a — /3| = 1. Then = 0.37, which 
certainly seems to be a low value. Accordingly, the states with a = and /3 = 1 seem 
very different. But what one can say about the states with a = 1000 and /3 = 1001? 
Their fidelity is low, but in many (of course, not all) cases replacing one of these two 
states with another hardly would cause a significant change in experimental results 
or their interpretation. On the other hand, the value of = 0.9 corresponds to the 
difference ja — /3| = 0.325, and it is not so obvious now, whether the coherent state 
with a = 0.325 can be considered as a good approximation of the vacuum state? These 
examples show that the fidelity alone not always can be used as an adequate quantity 
characterizing the degree of 'similarity' ('closeness') of quantum states. In many cases 
some other quantities should be taken into account in addition to the fidelity. 

One of the most important physical parameters is the energy. Therefore it seems 
interesting to answer the following question: how large the energy difference between 
two quantum states with the given fidelity can be? If it is small and the fidelity is 
high, then one may have more firm reasons to say that two states are 'close'. On the 
contrary, if the energy difference turns out to be big, then hardly two states can be 
considered as 'close', even if the fidelity is rather high. 

One can easily see that there is no unique answer to the question put above 
in the most general case. Indeed, let us consider two states, \^) and |x) = 
y/F\tp) + Vl ~ F\^±), where the state is orthogonal to 1-0). The fidelity between 
these states equals J" = KV'lx)!^ — F. Since the space of states \'tp±) is very big (I 
have in mind the continuous variable systems, when the Hilbert space of quantum 
states is infinite-dimensional), one can expect that there exist such states \ip±) that 
result in arbitrarily large difference of mean energies in the states \ip) and |x), without 
changing their fidelity (this is obvious if j-^) is the Fock state). 

But if one restricts somehow the space of admissible quantum states, then the 
above reasonings may become invalid (for example, there are no orthogonal coherent or 
squeezed states), and some universal (for the selected family) bounds can be found. It 
appears, in particular, that this is just the case for the family of Gaussian states. Since 
these states are frequently used in the contemporary theoretical and experimental 
quantum optics, and since the final results turn out remarkably simple, I believe that 
these results could be interesting for many readers. 



Relative energy difference of Gaussian states with a fixed fidelity 3 
2. Pure Gaussian states 

It is convenient to start from the simplest case of two pure squeezed quantum states 
described by the normahzed wave functions 

1 



{aj + ibj) {x — Xj) + ipjX 



where aj > 0, whereas real parameters bj, Xj and pj can assume arbitrary values 
(j = 1,2). For pure states, the fidelity ([T]) is reduced to the scalar product 
= |(V'i|'02)|^, which can be easily calculated: 

^ = 2v/aia2/Gexp(-C//G), (2) 

G= (ai+ 02)' + (61-62)', (3) 

U ={ai+ 02) {6pf+2 (0162 + 0261) {Sp){Sx)+[aia2 (ai + 02) + 0162 + aibj] {Sxf , (4) 

where Sx ~ X2 — xi and Sp = p2 — Pi- 

I suppose that the states {x) describe a quantum oscillator with unit mass and 
frequency. Then the mean energy of each state equals (assuming h = 1) 



-1 - [P] + ^3) /2 + (1 + a' + bj) I (4a,) . (5) 
The question is: how big can the difference — E2 — E\ be for the two states with 
different parameters but the fixed fidelity J-1 To simplify the formulas, it is convenient 
to introduce the following notation: a\ = a-2,~- a\ — aa, 61 — ac, 62 — 61 — a/3. Then 



6E = p{dp) + x{dx) + - [{Sp)^ + (Sx)^] + B, 

-c2)+/3(2c + /3)] 



B{a,c,a,l3) = 



a^ [a{l 



4a(l 



(7) 



Although the parameters a, /3, dx and Sp are limited for the fixed fidelity J^, the 
energy difference SE can assume arbitrarily big values, if the initial values x and p are 
very big. The same can happen if parameter a is very big or very small. Therefore it 
seems reasonable to study the limits of variations of the relative change of energy 

£ = 5E/Ei = E2/E1 - 1. (8) 

For the fixed values of 6p, 6x and Ei, the sum p{Sp) + x{Sx) is maximal if the vectors 
{6x,Sp) and {x,p) are parallel, and it is minimal for anti-parallel vectors. Therefore 
we have to find the maximal (positive) and minimal (negative) values of the functions 

±RS + B 



R^/2 + A 



(9) 



where A{a,c) = [l + (l + c^)] /(4a), R = ^x^ +p^, S = y/{Sxf + {5pf. 
Calculating the extremal values of ([9]) with respect to parameter R, we arrive at 
the problem of finding extremal values of functions 



under the constraint 
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\/(A + 2L)' + 8A ± (A + 2L) 



(10) 



4(1 + a) 



4(1 + a) + a2 + /32 



exp 



u{<p,a,c,a,l3) 
g{a,a,P) 



= T. 
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Here A = S^/A, L = B/A, g^a [(2 + + /S^] , 

u^{2 + a) sin^{ip) + [(1 + a)(2 + a + c^) + (c + pf] cos^{ip) 

+ 2a (2c + fi + ca) s'm{(p) cos{(p) 

and the new free parameter ip is the angle between the direction of vector (Sx, Sp) and 
the horizontal axis. 

2.1. The case of different displacements but identical squeezing parameters 

There are two special cases. The first one is a — (3 — (two squeezed states with 
different displacement parameters but identical variances of quadratures). Then 

£± = ±\[VA^ + 8A±a) , (11) 

, _ Sa'f[l + a^l + c^)r' 

sin^((/?) + (1 + c^) cos^((y3) + 2acsm{ip) cos{(p) ' 

where / = ln(l/J^). The maximal positive value of function ()lip and its minimal 
negative value (in the case of sign "minus") are achieved for the maximal possible 
value of the positive coefficient A. Looking for extrema of A as function of f one 
obtains Amax-ip = 4/k, where 

K=l + ^l~^, e(a,c) =4aVK(l + c2) +1]'. (12) 

Then equation ()lip results in the following limitations on the relative energy difference 
of two states with fixed parameters a and c (but arbitrary displacement parameters) 
for the given value of fidelity: 



<£ <^=^J- 1=. (13) 



One can notice an asymmetry between the left-hand and right-hand sides of inequality 
(|13p . Its origin is in the definition (jS]) , which is asymmetrical with respect to the states 
\tjji) and \4'2)- Due to this definition one has the restriction £ > —1 for negative values 
of £ (since the total energy of the state IV'2) is positive), whereas there are no bounds 
for positive values of £. The asymmetry in ([T3| disappears if one introduces the 
symmetric relative energy difference 

3^ = %=^ = ^. (14) 
Then both the inequalities in ([T3|) lead to the inequality 

y < V^- (15) 

One could normalize the difference |i?2 — £^i| in (ITi)) not by the 'geometrical' mean 
value ^/EiE2, but, say, by the usual 'arithmetical' mean value (Ei + E2) /2, but 
the resulting formulas are much more complicated and less attractive than the simple 
inequality (jl5l) . It would be interesting to know, whether there are some deep physical 
or geometrical reasons for choosing the geometrical mean values instead of arithmetical 
ones, or this is some mathematical artefact. 

In the special case of a = 1 and c = (when k = 1) we obtain the following 
exact relation between the maximal symmetric relative energy difference and the 
fidelity between two coherent states: 

3^1-") = V2ln{l/T). (16) 
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The inverse relation gives the maximal fidelity between two coherent states for the 
given value of y-. 

T^:^^^e^p{-yy2). (17) 

For arbitrary squeezed states with equal squeezing coefficients a and c, but 
different displacement parameters, one should take into account that coefHcient ^ 
goes to zero (and k — >■ 2) for highly squeezed states: either if a » 1 (strong coordinate 
squeezing), or a ^ 1 (strong momentum squeezing), or |c| — > cxd for a fixed a (strongly 
correlated states). For this family of states (labeled by the supescript S) we obtain 
the relations 



:y(f) = 2 VRT77), j-^l = exp (^y^A) . (is) 

2.2. Undisplaced squeezed states 

Another special case is (5 = 0. Then we have the constraint 

/il + a)^4D, D = {l ~ J^'') / (19) 

Formula ([9]) shows that the maximal relative energy difference £ in this case can be 
achieved for undisplaced squeezed states with R = 0. It can be written as 

£ = £c, + 2x{K + Mc), (20) 

where x = I [l + a'^ (l + (?)\ does not depend on a and /3, whereas £q,, K and M 
do not depend on a and c: 

£« = , K = 2D-£a, M ^ 



1 + a 1 + a 

For a fixed parameter c, the coefficient x varies from to (l + c^) ^ . Consequently, £ 
can vary between £a and £c — £a + 2{K + Mc) / (l + c^) . The extremal points of £c as 
function of variable c are given by the roots of the equation Mc^ + 2Kc — M = 0, i.e., 
Mc± = —K ± V + . These roots give us immediately the following extremal 
values of function £c{c): 



£c± =£»±( VK^ + AP±k) =2D± 2^/D{1 + D 



Surprisingly, these values do not depend on a or /3. They should be compared with 
the values of £a corresponding to maximal and minimal possible values of a for the 
fixed /3, namely a±(/3) — 2D ± ^y4D{l + D) — (3'^. Obviously, the minimal (negative) 
value of a_(/3) (achieved for /3 = 0) gives the maximal (positive) value of fa, whereas 
the maximal (positive) value of «+(/?) (also achieved for /? = 0) yields the minimal 
(negative) value £a. One can verify that these two extremal values coincide exactly 
with £c±. Thus we arrive at the inequalities 



<£<- , (21) 



y/l-T^ + l 1 - VI - -F^ ' 

which are equivalent to the following relations between the maximal symmetric relative 
energy difference ([T4| for the given fidelity and the maximal possible fidelity 
J^max for the given value of 3^: 



y -^-^HH J. 1 (22) 
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In the limit of Xji <C 1 or 1 — <C 1 we can write 

yrn ~ ^8(1- Tma. « 1 " yV&. (23) 

It is easy to see that 

-V!(coh) , y{S) ^ -V, -picoh) ^ -piS) , X- (24) 

•y-m ^ -'m ^ •^rm ■'max ^ ■'max ^ •'m.ax- V^^/ 

My conjecture is that the upper bounds ((22)) hold, as a matter of fact, in the most 
general case, when all three variations (a, /3 and d) can be different from zero. Many 
numerical tests made for different sets of parameters 14 confirm this conjecture: in 
all the cases the fidelity calculated by formula ^ appeared smaller than the maximal 
value given by with the parameter y calculated by means of equations (O and 
(|14p. But I did not succeed to find an analytical proof, except for the simplest (but, 
perhaps, the most important) special case when the difference 1 — J-" = e is small. 
Assuming that all variations (a, /?, 6x and Sp) are small and taking into account only 
the leading terms with respect to these variables, one can replace the exact formula 
for the relative energy difference £ by 

f = [pSp + xSx + ^ [a (1 - c") + 2c/3] - 

and to maximize this function under the constraint 

a^+/3^ + - \Sp^ + a^(l + c^) Sx^ + 2ac6pSx] = 8e. 
a 

This problem can be solved with the aid of the Lagrange multiplier, and the result 
coincides with (|23p . 



3. Mixed Gaussian states 



Since the extremal values of the fidelity and relative energy difference for pure 
Gaussian states are observed for zero displacement parameters [equation ([24])]. it 
seems reasonable to suppose that the same is true for mixed Gaussian states, as well. 
Therefore let us consider mixed homogeneous Gaussian states described by the density 
matrices (kernels of the statistical operators pk in the coordinate representation) 



Pkix,y) 



■ exp 



-i (ofe + ibk) - ^ (flfe - i^k) + akCkxy 



(25) 



where a^, bk and Cfe are real numbers obeying the inequalities Ok > and < Cfc < 1. 
The states ([25]) are normalized as follows: Tr (pk) = J Pk{x, x) dx — 1. The parameters 
^fc are responsible for the 'quantum purity': 



Aife = Tr {pI) 



'1-Cfc 



, 1 + Cfc- 

The mean energy in the state (pSjl equals 
l + aHl-C] 



Ck 



^1 



(26) 



(27) 



4a, (1-0) 

The calculation of the fidelity ([Ij between mixed Gaussian states is reduced to a chain 
of Gaussian integrals, since the square root of the Gaussian positive definite operator 
is also a Gaussian positive definite operator, whose kernel can be found using the 
scheme exposed, e.g., in [15] . Following this scheme one can arrive after some tedius 
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calculations to the following generalization of formula ^ for homogeneous mixed 
Gaussian states: 



^^^J^ , (28) 



^JG - (aiCi - a2C2)^ - 2Vaia2CiC2 



where the coefficient G is given by equation ([3]). Equivalent formulas for different 
parametrizations of the Gaussian states were found, e.g., in [5l fT5l - [T8] . 

3.1. Two .states with identical fixed purities 

Let us consider first the case of two states having the same fixed value of purity: 
Ci = C2 = C = const. Then one has to find the extremal values of the relative energy 
difference 5(a,/3) under the constraint 

"2(l_Sll£.,o. z,.<l^a-oii-c+^(i+oi^(29) 

The function £ in this case has the same form as in (PO]) . but with 

and unchanged coefficients £a and M. Using the same scheme as before, one can 
obtain the formula 



yra{FX) = yJ\-^[l + :F^a^-F)]. (30) 



Obviously, the right-hand side of ()30p decreases monotonously when C increases from 
to 1. Consequently, the most strong bounds on 3^,„ and Tmax take place for pure 
quantum states (C = 0), and they are given again by equation (|22l) . In particular, for 
C — >■ 1 ('supermixed' states with the purity /i — >■ 0) we have 



3;r"(-^)-2\/i-^, j---(3;) = ^-1^. (31) 



For small values of 1 — J-" and y 



3^r" « 2^/^~F, « 1 - :vV4. (32) 

3.2. Deviations from the pure state 

Another case which can be treated analytically is Ci = and C2 = C ^ 0. Then the 
constraint has the form 

(l + a)(l-C) ^^^^ 
The function £ can be written again in the form (|20|) with the same coefficient x^ but 
three other coefficients arc different: 

"~ (l + a)(l-C) "7^" (l + a)(l-C)' " (l + a)(l-C)" 

Then the same scheme as before leads to the formula 

y^n = (^^2(1-0 (i + c) --F2(i + c)' + y2(i-c) (i-c)-T^[i-cf^ 
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where C, = — C^. For C <C 1 we have 



3^™«- (^2^1--F2-C + ^(2J-2-l) + Cj . (34) 

The right-hand side of (|34p is obviously smaller than the pure state bound ((22|) . The 
maximal possible value of C, for the fixed T follows from equation ([55]) with a ~ /3 = 0: 

Cma. = 2^-2 (l - VI --^'(l--^')) . If e = « 1, then Cma. = 2£-5eV2 + - • •. 
Putting this critical value in p4p one can obtain the value ymiCmax) — 2e, which is 
much smaller than the pure state boundary ^/Se. 



4. Discussion 



The results obtained are illustrated in two tables below. The first one gives the 
maximal possible fidelities between two coherent states [equation P^ ]. two displaced 
pure squeezed states with arbitrary fixed variances [equation (|T8l) ]. two 'supermixed' 
homogeneous squeezed states with C ~^ 1 [equation (|3T|) ] and two arbitrary (pure or 
mixed) squeezed states with zero displacements [equation (22)] for fixed values of the 
energy ratio E2/E1. The second table gives the values of the maximal symmetrical 
energy difference and the corresponding maximal possible ratio {E2/ Ei)^_^ for the 
given value of fidelity between the most general Gaussian states [equation (l22t ]. 



E2/EI 


y 


•J inax 


■J max 


■psmix 

^max 


max 


F 


ym 


{E2/E,),^ 


3 


1.155 


0.51 


0.72 


0.75 


0.87 


0.999 


0.09 


1.09 


2 


0.707 


0.78 


0.88 


0.89 


0.94 


0.99 


0.28 


1.32 


1.5 


0.408 


0.92 


0.96 


0.96 


0.98 


0.95 


0.66 


1.91 


1.1 


0.095 


0.9955 


0.9977 


0.9977 


0.998866 


0.9 


0.97 


2.55 



It is worth emphasizing that there exist quantum states with the values of y or 
T which are arbitrarily close to the maximal ones found above (for the fixed value of 
the other parameter). In particular, J-max = \/8/9 ~ 0.94 for E2 = 2Ei or Ei = 2E2. 
This means that even 94% of fidelity cannot guarantee that two squeezed states are 
really close in energy to each other {in the absence of any additional information about 
the states). Explicit examples are two states with bj = xj = pj = and 02 = 2ai 
(then a = 1) or a2 = ai/2 (then a = -1/2), with the fidelity 94%. Taking ai = 1 (the 
oscillator ground state), one obtains the mean energies (£^1) = 1/2 and {E2) ~ 5/8, 
which are not too different (3^ = l/y/20 « 0.22). However, for the highly squeezed 
states with ai = 0.1 and 02 = 0.05 one has (Ei) — 2.525 and {E2) — 5.0125; hardly 
these two states with y « 0.7 can be considered as 'close'. On the other hand, if one 
knows, for example, that the first state is coherent with xi = 100 and pi = 0, and 
the second state is also coherent with p2 = 0, then T = 0.94 means that X2 = 100.35 
or X2 — 99, 65. In this case y ~ 0.007, which can be accepted as a small quantity, 
meaning that these two concrete states are 'close'. But to guarantee that the relative 
energy difference is below, say, 10% {E2 < 1.1-Ei) in the absence of any additional 
information, the fidelity must be higher than 440/441 w 0.998866. 

Formulas derived in this paper can be rewritten in terms of the Bures-Uhlmann 

distance \ll B = (2 — 2V^^ ■ But the corresponding expressions seem less 
attractive. For example, instead of (1^^ one obtains the following formula for the 



Relative energy difference of Gaussian states with a fixed fidelity 



9 



minimal Bures distance between two Gaussian states with the fixed relative energy 
difference: 



The mutual bounds on the fidelity and the relative energy difference can be 
derived also for other interesting families of quantum states. In some cases (e.g., 
for binomial and negative binomial states) analytical formulas can be found [19] . but 
in the most of cases (e.g., for superpositions of coherent states) this can be done only 
numerically. Perhaps, it would be interesting to obtain different bounds taking into 
account other parameters (besides the energy) characterizing quantum states, such as 
the degree of squeezing, Mandel's parameter, or something else. 
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